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We study the long term evolution of magnetic fields generated by an initially unmagnetized collisionless relativistic 
e~^e~ shock. Our 2D particle-in-cell numerical simulations show that downstream of such a Weibel-mediated shock, 
particle distributions are approximately isotropic, relativistic Maxwellians, and the magnetic turbulence is highly 
intermittent spatially, nonpropagating, and decaying. Using linear kinetic theory, we find a simple analytic form for 
these damping rates. Our theory predicts that overall magnetic energy decays like (u)pt)~'^ with g ~ 1, which compares 
favorably with simulations, but predicts overly rapid damping of short wavelength modes. Magnetic trapping of 
particles within the magnetic structures may be the origin of this discrepancy. We conclude that initially unmagnetized 
relativistic shocks in electron-positron plasmas are unable to form persistent downstream magnetic fields. These results 
put interesting constraints on synchrotron models for the prompt and afterglow emission from GRBs. 
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1. Introduction 

The prompt emission and afterglows of gamma-ray 
bm'sts (GRBs) may be manifestations of ultrarela- 
tivistic shock waves. These shock waves may be me- 
diated via the relativistic form of Weibel instability 
(Weibel 1959; Yoon and Davidson 1987; Medvedev 
and Loeb 1999; Gruzinov and Waxman 1999). The 
free energy from strong plasma anisotropy in the 
shock transition layer generates strong magnetic 
fields (with strengths comparable to the available 
free energy). However, these fields have very small 
spatial scales, i.e., the order of the plasma skin depth, 
c/ujp, where uop is the plasma frequency. These ini- 
tially small-scale B-ficlds must survive for tens of 
thousands to millions of inverse plasma periods to 
serve as this source of the magnetization for syn- 
chrotron models of burst emission and afterglows 
(Gruzinov & Waxman 1999; Piran 2005ab; Katz, 
Keshet, & Waxman 2007). Whether or not these field 
can is an open question. 

Numerical and analytic studies (Kazimura et al. 
1998; Silva et al. (2003); Frederiksen et al. 2004; 
Medvedev et al. 2005; Hededal et al. 2005; Nishikawa 
et al. 2003, 2005; Spitkovsky this proceedings) have 
elucidated the basic physics. The instability initially 
forms filaments of electric current and B fields, 
which then merge to inverse cascade magnetic en- 
ergy to larger scales, but only in the foreshock re- 
gion. When the B-fields reach the magnetic trap- 
ping limit (Davidson et al. 1972; Kato 2005; also see 
Milosavljevic, Nakar, & Spitkovsky 2006; Milosavlje- 



vic & Nakar 2006a), particle orbits become chaotic, 
disorganizing the filaments. The disorganized mag- 
netic fluctuations scatter their supporting particles, 
which isotropizes and thermalizes the flow, within 
tens to hundreds of skin depths (Spitkovsky 2005). 
The magnetic energy peaks in this layer at ~ 10-20% 
of the bulk plasma flow energy. 

However, present simulations have not deeply 
followed the flow into the downstream region to ex- 
plore the long term behavior of these B-ficlds. Thus, 
the question of the structure and long-term survival 
of the B-fields remains open (see for instance, Gruzi- 
nov & Waxman 1999; Gruzinov 2001ab; Medvedev 
et al. 2005). 

In this proceeding, we discuss recent work 
(Chang, Spitkovsky, and Arons 2008; hereafter 
CSA08) which shows that this magnetic energy must 
rapidly decay in the downstream medium. We first 
describe the basic features of the downstream plasma 
from our numerical simulations. We then calculate 
the evolution (decay) of the downstream plasma us- 
ing Vlasov linear response theory and then compare 
this evolution with simulations. While linear theory 
does reasonably well in estimating the decay rate of 
the total magnetic energy, it overestimates the damp- 
ing rate of shorter wavelength modes. We discuss this 
discrepancy as a result of magnetic trapping. Finally, 
we summarize our results. 
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2. Simulation Results 

Spitkovsky (2005, this proceedings) and Spitkovsky 
and Arons (in prep) describe a series of 2D and 3D 
simulations of relativistic shock waves in e^e~ plas- 
mas. These are Particle-in-Cell (PIC) simulations, 
using the code TRISTAN-MP. We simulate shocks 
by injecting cold relativistic plasma particles at one 
end of a large domain that reflect off a fixed conduct- 
ing wall at the other end. We use 2D boxes as large 
as 50,000 X 2048 cells with up to 1.35 x 10^° particles 
to study these shocks. We refer the reader to CSA08 
for additional details. 
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Fig. 1. Snapshot of a region from a large 2D relativistic shock 
simulation in the downstream frame, a) Density structure, nor- 
malized in upstream units, in the simulation plane, b) Mag- 
netic energy, normalized in terms of upstream energy of the 
incoming flow; eg = / Air^yinimc^ . A power law scaling, 
e^J'^ , was applied to stretch the color table to show weak field 
regions and is reflected in the colorbar. c) Plasma density aver- 
aged in the transverse direction as a function of the distance 
along the flow, d) Magnetic energy density averaged in the 
transverse direction, as a function of distance along the flow. 



downstream (Fig. Ic), which is the expected com- 
pression factor (Gallant et al. 1992; Spitkovsky and 
Arons, in prep). Current filaments show up as an en- 
hancement in the plasma density and magnetic en- 
ergy density in the foreshock (Fig. lab). The scale of 
the filaments grows towards the shock through merg- 
ing. 

At the shock transition layer, the filaments disor- 
ganize and become clumps of magnetic energy. These 
magnetic clumps lose intensity the further down- 
stream they are from the shock (Figure lb). We also 
find that the particle distribution function changes to 
an isotropic (in the downstream rest frame) thermal 
population, i.e., the difference between the perpen- 
dicular and parallel momentum is extremely small, 
< 1%. As Figures 1 b and d suggest, the magnetic 
fields decay in the downstream region of the shock. 



3. Downstream Evolution of Magnetic 
Turbulence 

We now attempt to analytically understand this de- 
cay of magnetic turbulence. The simulations show 
that the downstream plasma is isotropic and the 
downstream particle distribution function is well de- 
scribed by a relativistic Maxwellian. We will also as- 
sume the downstream field amplitudes are so small 
that particle orbits are almost straight lines. 

We begin by deriving the linear plasma response 
is determined by the plasma susceptibility, x (Stix 
1992). We evaluate the susceptibilty for distribution 
functions that are isotropic in two and three dimen- 
sions. Details can be found in CSA08. We set = 0, 
because of the non-propagating nature of the mag- 
netic clumps, which we infer from the simulations. In 
the long wavelength limit, i.e., k <^ i-Op/c, where k is 
the wavenumber, we find: 



Figure 1 shows the snapshots of density and 
magnetic energy from a typical 2D simulation. Co- 
ordinates are in units of the upstream skin depth, 
c/ujpi. In the simulation shown, the upstream flow 
moves to the left with 71 — 15. 

Our simulations are large enough to permit the 
complete development of the shock and show the 
main features of contemporary coUisionless shock 
simulations. For instance, we see the factor of « 3.13 
increase in density between the upstream and the 
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Note that the 2D and 3D results only vary by a 
numerical factor. Long wavelength modes have the 
same qualitative behavior in two and three dimen- 
sions. 

The plasma susceptibility (eq. [1] ) can be utilized 
to calculate the evolution of an initial field of fluctu- 
ations. We refer the interested reader to CSA08 for 
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details, but the result is 

d\6Bk\^ ^ 
dt 



(2) 



where 7fe = (fcc)^tj (Attx) ■ The asymptotic 
forms of X from equation (1) for 2D and 3D gives 

2D 

7fc = <; .7..,3 _ . (3) 



4 |A.;c|' 



3D 



Note the strong cubic k dependence on the decay 
rate. Short wavelength modes rapidly damp, but 
longer wavelength modes can persist. 

We now compare these expectations to the nu- 
merical simulations. We take the Fourier transform of 
SB from our 2D numerical simulations from a down- 
stream region behind the shock front at a; = a;o. We 
evolve these spectra for 450 (red), 900 (green), and 
1350 LUp^ (blue) using equation (3). We compare this 
analytically evolved spectra to Fourier transformed 
snapshots taken from our numerical simulations at 
these times in Figure 2, where xq = 8A0c/ujp. While 
theory and simulation agree at very low wavenum- 
ber (kyc/ujp < 0.2), theory overpredicts the cutoff in 
power at larger k. The discrepancy may be due to 
magnetic trapping (see §4). 
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Fig. 2. Spectral evolution of magnetic field from the slice at 
SiOc/up. Initial field spectrum (black solid line) is plotted after 
450 ojp^ (red), 900ajj7^ (green), and 1350 ojjj"^ (blue) based on 
simulation data. Dashed curves represent analytic evolution 
of the initial field. 



Since total B-field energy is dominated by long 
wavelength modes, we use equation (2) to find a sim- 



ple decay law for the total B-field. Again we refer the 
interested reader to CSA08 for details, but if the ini- 
tial spatial spectrum is a power law in wave number 
|(5i?fep oc /c^P, then the B-field should decay like 



5B^ oc i-2(p+i)/3. 



(4) 



For a shock moving at constant velocity, we have 
a^poak-a; oc t. Hence SB"^ oc (a^pcak — a;)"^^^"*"^^/^. Our 
numerical simulations are extremely suggestive that 
the magnetic energy density follows the a SB'^/Stt oc 
^-2/3 decay expected for an initially flat magnetic 
spectrum at early times {p = 0), then steepening to 
a t^^ decay at later times {p — 1/2) as shown in 
Figure 3. We have analyzed additional simulations 
with a large transverse spatial scale and they sug- 
gest p — 0. This difference in the index of the decay 
law expected from the theory and measured from the 
simulations may also be due to magnetic trapping 
(see §4; see Gruzinov 2001b for an alternate expla- 
nation and CSA08 for a rebuttal). 




Fig. 3. Magnetic energy density (in units of upstream kinetic 
energy) as a function of position downstream of the shock. A 
broken power law proportional (x — ajpeak)"^''^ fits well at 
early times, but a (a; — ipeak)"'^ power law fits better at later 
times. 



4. Magnetic Trapping 

While our simulations and theory are consistent with 
one another in suggesting overall power law decay t~'' 
with q ^ 1, they disagree for short wavelengths. This 
discrepancy may arise from the nonlinear effects of 
magnetic trapping. Our simulations also show that 
not all particles follow straight line trajectories that 
are weakly perturbed, but some are partially trapped 
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and strongly deflected. By following test particle or- 
bits in our simulations, we find that the Larmor radii 
of many of the test particles are of the same order 
of the sizes of these clumps or smaller. These strong 
departures from weakly perturbed particle dynamics 
may be the cause of the decreased damping at large 
wavenumber found in the simulations. It may also 
modify the overall decay away from t~'^/^ decay law 
expected from simple linear theory. 

5. Discussion 

We have studied the downstream evolution of mag- 
netic turbulence in the context of a collisionless e+e~ 
shock both analytically and numerically. Our simu- 
lations show that the downstream region consist of 
nonpropagating magnetic clumps embedded in quasi- 
homogenous medium where the background parti- 
cle distribution function is an isotropic Maxwellian. 
In such a background, we showed that magnetic en- 
ergy will decay like t~'^ with g ~ 1. However, linear 
theory overpredicts the decay rates at short wave- 
lengths compared to simulations. Magnetic trap- 
ping may play an role in resolving this discrepancy. 
Rapid field decay puts severe constraints on GRB 
emission mechanism, but they may not be incon- 
sistent with GRB observations (see Pe'er & Zhang 
2006). Finally, if ion-electron collisionless shocks 
reach roughly equipartition with each other as sug- 
gest by recent large scale simulations (Spitkovsky 
2008), they would reproduce the physics of the 
shock and their B-flelds would decay as well. 
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